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ABSTRACT 

Microlensing surveys, which have discovered about a dozen extrasolar planets to 
date, have focused on the small minority of high-magnification lensing events, which 
have a high sensitivity to planet detection. In contrast, second-generation experi- 
ments, of the type that has recently begun, monitor continuously also the majority 
of low-magnification events. We carry out a realistic numerical simulation of such 
experiments. We simulate scaled, solar-like, eight-planet systems, studying a variety 
of physical parameters (planet frequency, scaling of the snowline with stellar mass 
-Rsnow ot Af), and folding in the various observational parameters (cadence, exper- 
iment duration), with sampling sequences and photometric error distributions taken 
from the real ongoing experiment. We quantify the dependence of detected planet yield 
on cadence and experiment duration, e.g., the yield is doubled when going from 3-hour 
to 15-minute baseline cadences, or from an 80-day-long to a 150-day-long experiment. 
There is a degeneracy between the snowline scaling index s and the abundance of 
planetary systems that can be inferred from the experiment, in the context of our 
scaled solar-analog model. After 4 years, the ongoing second-generation experiment 
will discover of the order of 50 planets, and thus will be able to determine the fre- 
quency of snowline planet occurrence to a precision of 10-30%, assuming the fraction 
of stars hosting such planets is between 1/3 and 1/10, and a snowline index in the 
range s = 0.5 to 2. If most planetary systems are solar analogs, over 65% of the de- 
tected planets will be 'Jupiters', five in six of the detected anomalies will be due to a 
single planet, and one in six will reveal two planets in a single lensing event. 

Key words: surveys - gravitational lensing: micro - methods; statistical - binaries: 
general - planetary systems - Galaxy: stellar content 



1 INTRODUCTION 

Among the techniques currently used to discover extrasolar 
plan ets, microlensing has some unique capabilities (see, e.g. 
I Goul d ot al. 2010). Apart from the fact that microlensing 
probes for planets around stars over a large volume of the 
Galaxy, among both the bulge and the disk population, and 
can also discover free-floating planets (|Sumi et aLlboill ). it 
is sensitive to planets at projected host separations that 
are distinct from those probed by other methods, namely 
about 1 to 10 AU. This corresponds to the area beyond the 
"snowline" where gas and ice giants are likely to form. Mi- 
crolensing is thus complementary to other methods, such 
as radial velocity or transit surveys, which are most sensi- 
tive to planets in close orbits around their host stars, and 
to direct imaging of planets, which probes larger separa- 



tions (e.g. ICrepp fc JohnsonI 120111 ). To date, however, only 
abo ut a dozen p l anets have been discovere d via microlens- 
ing JBond et all |2004|. lUdalski et al.1 1 2005 |. iBeaulieu era! 
Gaudi et al. 2008. Bennett et al 
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iBatista et al.ll2011^ ■ out of the to- 



tal ~550 known planeto and the ~1200 candidates rec ently 
announced by the Kepler mission (jBorucki et al.ll201ll ). 

The paucity of microlensing-discovered planets is a con- 
sequence of the surveying strategy, largely followed to date, 
which has focuse d on high-magnific ation lensing events. 
As summarized in lGould et al.l ((20101), all currently known 
microlensing-discovered planets were first discovered by one 
or both of two wide-fi eld, low-cad ence microlensing surveys, 
OGLE (|Udalskill200"9l ) and MOA (|Sako et al.ll2008l ). Events 
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whose initial light curves suggested high peak magnifica- 
tions triggered global networks of telescopes for photomet- 
ric (and sometimes spectroscopic) follow-up of these specific 
events. Indeed, among the 12 known microlensing planets 
(11 systems, as OGLE-06-109L has two planets) eight were 
discovered in events with amplifications greater than 100, 
and three had medium magnifications (greater than 10). 
All of these events had long event time scales (40-60 days) , 
whereas the typical event time scale is ~20 days (see Section 
[S] below), and high cadences we re therefore not esse ntial. 
Only one event, OGLE-05-169L (|Beaulieu et al.l 120061 '). was 
a low-magnification event, with an amplification of 3.0 and 
a relatively short event time scale of 11 days. The source, 
however, was a bright star with a baseline magnitude of 
I = 14.3, which permitted triggering of the follow-up net- 
work nonetheless. All of the above events had peak /-band 
magnitudes brighter than 16.7. 

Given the limited observational capabilities to date, the 
high-magnification follow-up strategy has been the optimal 
choice. Bright events can be followed up with high cadence 
using small telescopes, and these events are fortuitously also 
those with the highest sensitivity to planets. However, only 
1% of all microlensing events have magnifications > 100. 
Since the surveys of the Galactic bulge have discovered 
roughly 700 events per year in recent years, and apparently 
only a fraction of all stars host planets near their snowlines, 
only a couple of planets per year have been discovered, re- 
sulting in the current microlensing planet tally. 

Microlensing survey statistics have been used in the 
past to constrain the frequency and p r operties of snowline - 
regi on planets bylAlbrow et alj (|200ll ') , iGaudi et aJ] (|2002l) , 
and lTsapras et al.l (|2003t ). More recentlv. [Sumi et all (|2010l ') 
derived the mass-ratio distributions of cold exoplanets, and 
found that Neptune-mass planets are th ree times more com - 
mon than Jupiters beyond the snowline. iGould et al.l (|2010l ) 
find from the high-magnification events that the frequency 
of Solar-like systems around lens stars is about 1/6. How- 
ever, two main problems complicate these estimates. One 
is the small-number statistics they are based on. The other 
is the complex social process that determines which events 
are followed up, and how, a process which then needs to 
be modeled in order to estimate the intrinsic planet proper- 
ties. iGould et al.l (|2010l ) argue that the very chaotic nature 
of the process results, in the end, in a randomization that 
minimizes selection effects and gives a high level of com- 
pleteness. Nevertheless, a more controlled microlensing ex- 
periment would provide confidence that this is indeed the 
case. 

It has been clear for some time that, to make progress 
and overcome these problems, a "second gener ation" (here- 
after generation-II) microlensing survey (e.g., iGaudi et al.l 
[2OO9I) must be initiated, one in which a large fraction of 
all ongoing microlensing events toward the Galactic bulge - 
both low and high magnification - are monitored continu- 
ously with a cadence that is high enough to detect planetary 
perturbations in the light curves. Although only a fraction 
of low-magnification events will reveal planets, even if they 
are present around the lens star, the numbers of such events 
are much larger than high-magnification events. In the bal- 
ance, such a survey, over several years, could yield a sample 
of planets that is competitive with those from other tech- 
niques. The requirements are straightforward - a network of 



Im-class telescopes situated so as to allow 24-hour coverage 
of the bulge, equipped with degree-scale imagers. In June- 
July 2010, we carried out a 6-week pilot version of such a 
generation-II experiment. Starting April 2011, we have be- 
gun a full generation-II experiment. In this paper, we per- 
form numerical simulations of the possible outcomes of this 
and similar experiments. 

We begin our study by considering a variety of physi- 
cal possibilities for the abundance of planetary systems and 
their properties. We then ray trace through the systems to 
generate a library of microlensing light curves with the cor- 
rect distribution of observed properties, such as peak mag- 
nification, peak magnitude, and timescale. The simulated 
light curves are then "observed" using the real observational 
parameters of our ongoing experiment, including sampling, 
weather gaps and photometric noise. Finally, the data are 
processed through a realistic detection pipeline, to recover 
the detection efficiency and the expected number of detec- 
tions for each combination of physical and observational pa- 
rameters. The results we present here can therefore serve 
to predict the power of an experiment to distinguish among 
models, and to optimize the experiment's yield. We note 
that, although generation-II experiments will discover many 
well-characterized planetary events, our focus is on the mere 
detection of planetary anomalies, even when they may end 
up not being unambiguously modeled. 

Our methodology builds upon, and extends previous 
studies of the pl anet detection efficienc i es of microlens- 
ing surveys (e.g., ^ Gaudi fc Sacketd l200d: Bennett fc Rhig 
2002; Bennett 20041: IYoo et al.l I2OO4I : iDominik et al.l I2OIOI : 
Pennv et al.ll2011l ). However, to our knowledge, ours is the 
first study to combine the following important elements in 
such a simulation: we consider full 8-planet solar-analog sys- 
tems; we investigate the effect of the snowline radius of those 
systems, for a range of scaling dependences on mass; for our 
simulated light curve generation, we use the observed distri- 
butions of lensing event timescales, impact parameters, and 
peak magnitudes, and the distributions of lens mass, dis- 
tance, and velocity implied by Galactic models; we use the 
real sampling patterns and photometric error distributions 
(including outliers) of the ongoing generation-II experiment 
that we simulate. Eventually, we will use the tools we de- 
velop here to constrain the actual physical parameter space 
of exoplanets probed by microlensing. 

In Section O we present the observational setup of our 
generation-II network. In Section [3] we describe the simula- 
tions, starting with the light curve generation process, and 
up to the detection procedure. In Section [3] we discuss the 
results we obtain as a function of the various input param- 
eters. A summary of our results is given in Section [S] 



2 THE GENERATION-II NETWORK 

The generation-II microlensing experiment is, in essence, a 
continuous-monitoring, high-cadence survey with a field of 
view of 8 deg^ centred on the Galactic bulge. The experiment 
network currently includes four telescopes: 

(i) The Im-telescope at Wise Observatory in Israel 
with the LAIWO ca mera, with a 1 deg^ field of view 
(jGorbikov et al.l I2OI0I ) . The cadence for each of the eight 
Wise fields is ~30 minutes. 
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Table 1. The generation-II network 



Group (site) 



Aperture FOV 

[m] [deg2] 



Cadence 

fminl 



High cadence 
area [deg^] 



Wise (Israel) 


1.0 


1.0 


30 


8.0 


OGLE (Chile) 


1.3 


1.4 


15 (45) 


4.2 (7.0) 


MOA (NZ) 


1.8 


2.2 


15 (45) 


13.2 (13.2) 


PTF (California) 


1.2 


7.8 


40 


7.8 



Note — numbers in parentheses are the cadences and areas of the 
lower-cadence OGLE and MOA fields. 




Figure 1. Gcncration-II network fields: red — Wise, blue 
OGLE-IV, green - MOA-II, magenta - PTF. 



(ii) The Las Campanas Observatory 1.3m Warsaw Uni- 
versity telescope in Chile with the OGLE-IV camera, with a 
1.4 deg^ field of view (Udalski 200a' ). The OGLE-IV phase 
has three extremely high-cadence fields observed every ~15 
minutes and five high-cadence fields observed once per ~45 
minutes. 

(in) The 1.8m MOA-II telescope in New Zeala nd with the 
M0A -cam3 camera, with a 2.2 deg^ field of view IjSako et al.l 
120081 ). The MOA-II phase has six extremely high-cadence 
fields (every ~15 minutes) and six high-cadence fields (every 
~45 minutes) ; most of these latter fields are outside the full 
network footprint. 

(iv) The Palomar Observatory 1.2m Oschin telescope in 
California with th e PTF camera, with a 7.8 deg^ field of view 
IjLaw et al.l 12003 ) and a cadence of ~40 minutes. This tele- 
scope participates only during one month every year, when 
the network has the highest continuous coverage. Further- 
more, the PTF footprint overlaps with only 4 deg'^ of the 
Wise fields. 

The network parameters are summarized in Table [l] Figure 
[l] shows the region of the Galactic bulge covered by the 
network. The 8 deg^ covered by Wise, OGLE, and MOA 
include the fields with the highest microlensing event rates 
in previous OGLE and MOA seasons. 



3 EXPERIMENT SIMULATION 

In this section, we describe the various stages of our nu- 
merical simulation of the generation-II microlensing planet 
survey. 



3.1 Source- lens configuration 

A point-mass microlensing event is characterized by its an- 
gular Einstein radius, 9e'- 



AGMl Ds - Dl 



(1) 



c2 DsDl 

where Ml is the mass of the lens star, Ds and Dl are the 
distances to the source and lens stars, respectively, G is the 
gravitational constant, and c is the speed of light. The time 
scale of the event, Ie, is the time it takes the source star to 
cross the angular Einstein radius due to the star's relative 
transverse velocity, vt- The angular position of the source 
star at time t, neglecting parallax (due to Earth's acceler- 
ation), orbital motion of the lens due to a companion, and 
"xallarap" orbital motion of a binary source is: 



u(t) 



< + 



t-to 

tE 



(2) 



where u{t) is the angular separation between the source and 
the lens in units of Oe, wo is the impact parameter - the 
minimum angular separation between the source and the 
lens, and to is the time at minimum separation. Among ts, 
«o and to, ts contains the interesting physical properties 
of the source-lens configuration (distances, relative veloc- 
ity and lens m ass). Our first s tep is to model those physi- 
cal properties. iDominikI (120061 ) has used Galactic structure 
models to find the probability functions for the lens mass, 
lens distance, and transverse velocity, as a function of the 
event time-scale, tE, in the direction of the Galactic bulge 
through Baade's window (Galactic longitude and latitude 
[l,b] — [1°,— 3.9°]), for a source star at Ds = 8.5 kpc (see 
his figure 5). We create our simulated sample of lenses by 
choosing ts from the observed distribution of all OGLE-III 
(2002-2009) events (Fig. [I}, as listed in the OGLE databasfE 
and, for each tE, drawing Ml, Dl and vt from the iDominik 
(2006) probability functions. 

The impact parameter, no, should, in principle, have a 
uniform distribution. In practice, however, due to the fiux 
limits of the actual observations, there is a selection effect 
that, for intrinsically faint sources, favors the detection of 
high- magnification events (i.e., small uo). In line with our 
objective of making our simulations as realistic as possible, 
we therefore choose also the impact parameter from the dis- 
tribution of all OGLE-III events (Fig. [2]). We confirm that, 
for OGLE events where the source was unblended, implying 
it was always above the limiting magnitude, the distribution 
of Wo is consistent with being uniform (up to Poisson fluctu- 
ations). With this, we have the full source-lens configuration 
of an event. 
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Figure 2. Left: Observed distribution of event time scales, Ie, for all OGLE-III (2002 — 2009) events with uq ^ Ij obtained from the 
OGLE database. The solid curve is a lognormal fit, A^ oc exp(— ^ °^\^ 'i '^ ), with /x = 1.28 and a = 0.43. Right: Distribution of impact 
parameter, uo, for all OGLE-III events with no ^ 1 (light-gray). Due to selection favoring the detection of highly magnified faint events 
over weakly magnified faint events, there are more events at smaller impact parameters. The dark-gray histogram shows the distribution 
only for unblended (i.e. relatively bright) events, which are seem to be distributed uniformly, as expected. 



3.2 Planetary systems 

The characteristics and diversity of "normal" extrasolar 
planetary systems are still virtually unknown. Although 
there is an infinite variety of ways in which planetary sys- 
tems could, in principle, be structured, with planets of all 
masses at all radii, real multiple planetary systems obey 
physical constraints, resulting from formation history and 
stability criteria. A Copernican approach is therefore to 
assume as a working hypothesis (to be tested by observa- 
tions) that the Solar system is a typical planetary system, 
as hinted by the microlensing discovery of a scaled-down 
Sun- Jupiter-Saturn syste m in OGLE-06-109L (JGaudi et al.l 
120081 : iBennett et al.ll201Cl ). Indeed, a major advantage of mi- 
crolensing as a planet-hunting technique is that it is the most 
sensitive method available for the detection of the snowline 
region of such analogs (Gould et al.ll2006l . lGould 2009 ). We 
thus choose our simulated planetary systems to mimic the 
Solar system, with eight coplanar planets in circular orbits, 
having the same mass and distance ratios as the Solar sys- 
te m. The ste llar mass, Af_L, as noted above, is drawn from 
the iDominikI (J2006) probability function. 

The location of the snowline at the time of planet for- 
mation is thought to set the border between rocky planets 
and gas/ice giants, and therefore the dimensions of a plane- 
tary system may scale with the snowline radius. The snow- 
line radius itself will depend on the stellar mass, in a way 
that can be parametrized as a power-law, fisnow oc M" . The 
value of the index s has been debated (e.g., iLecar et al.ll2006l . 
[Kennedy et~ai]|2006l . Ikennedv fc Kenvonll2008l ). For exam- 
ple, considering just the equilibrium temperature around 
main-sequence stars, s ~ 2 would be expected. However, var- 
ious effects at the time of planet formation in a protoplane- 
tary disk around a pre-main-sequence star could modify this 
value. In our simulation we explore the range 0.5 ^ s ^ 2. 
The zeropoint for the scaling is Jupiter's orbit in the solar 
system (5.2 AU). 



Finally, for each simulated event, the inclination to our 
line of sight of the planetary system, the phase of each planet 
in its orbit, and the source trajectory angle with respect to 
the orientation of the orbital plane of lens system on the sky, 
are chosen at random. 



3.3 Light-curve generation 

For a lens that is a star with a planetary system that in- 
cludes n — 1 planets, the lens equation that relates the source 
position, p , and the image positions, 9 , in the lens plane 

is: 



9%,L 



E' 



(3) 



where 8e,l is the angular Einstein radius of the lens star, 
qi is the mass ratio between the zth planet and the lens star 
{qi — 1 for the lens star) and 9i is the position of the ith 
planet. The magnification is the ratio of the summed solid 
angles subtended by all the images to that of the solid angle 
subtended by the source. 

In principle, solving the lens equation for the image po- 
sitions, given a source position, requires solving a polyno- 
mial of degree n^ + 1 although some of those solutions are 
not physical, and hence the number of variables can be re- 
duced to 2n (|Wittlll990r) . Given the large number of planets 
per star that we wish to investigate in our simulations, we 
choose to generate lense d images via inverse ray-shooting 
(|Schneider fc Weisg|l988l ) and thus solve for the magnifica- 
tion A{t) of the source flux. In this brute-force approach, 
we divide the lens plane onto a grid and use the lens equa- 
tion directly to map the lens plane onto the source plane. 
In order to achieve fast performance and high accuracy, we 
use an adaptive grid that increases the lens plane resolution 
at the positions of the images. We can thus also account 
for finite source effects (|Nemiroff fc Wickramasinghd 119941 ) 
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and for the limb-darkened profile of the source star. For sim- 
plicity, in this work we use a constant source size of IRq, 
and a uniform source surface brightness. Giant sources with 
radii of ~ 5-10 Rq, which can be a substantial fraction of 
the sources, can eliminate small-deviation detections due to 
their larger area but, on the other hand, have a larger effec- 
tive cross-section for planet lensing. Assuming a rough bor- 
der in unlensed source magnitude of /= 17.5 mag between 
giants and dwarfs in the bulge we find that about 20% of cur- 
rent OGLE microlensing events have giant sources. To gauge 
their effect, we found that, even if all the sources in our sim- 
ulations, described below, are assumed to be giants, fewer 
than 5% of detections are missed. Figure [3] shows several 
examples of planetary configurations (left panels) and the 
microlensing light curves that they generate (central pan- 
els). 

3.4 Event sampling 

The last step in our procedure is to simulate realistic ob- 
servations of the light curves that we have generated. We 
include the following effects. First, we sample each event in 
a way that mimics the real observations, including gaps due 
to weather, full moon, and technical problems. Specifically, 
we use the actual sampling patterns of the generation-II net- 
work from the pilot campaign. The cadence of the sampling 
is crucial for planet detection because the time scale for 
the duration of the planetary devia tion is relatively short, 
^planet ~ y^-^B l|Gould fc Loeblll992h . A baseline magnitude 
needs to be assigned to the event, which may be below the 
detection limit of a given telescope. Frequently, the source 
is blended with other stars close to the line of sight. Pho- 
tometric errors added to each measurement need to include 
Poisson noise from both the source and the stellar back- 
grounds. In a real experiment, however, difference-image 
analysis (DIA) is usually performed in order to achieve accu- 
rate photometry in crowded fields. The photometric errors 
on the difference residuals are often under-estimated and 
there may be residuals due to imperfect subtractions, re- 
sulting in outlier points in the light curves, and these errors 
need to be simulated as well. 

We include all of these effects as follows. Every simu- 
lated event is matched to a real event, drawn from the list 
of all past OGLE-III events. For each such event, as noted 
in Section l3.ll we use the actual Ie, and also the source 
base magnitude and blending fraction. Each light curve is 
sampled with a pattern of "minimal" cadences (from 15 min- 
utes up to 3 hours) that are later corrupted by the actual 
sampling patterns of the 2010 pilot experiment, including 
measurements lost to weather, full moon, and so on. Fur- 
thermore, our sampling reflects the different sampling pat- 
terns of each observatory. Thus, for example, changing the 
minimal cadence of a simulated run from 15 minutes to 30 
minutes does not affect the PTF sampling pattern, which 
already has a 40-minute cadence to begin with. The total 
duration of the experiment per bulge observing season is 
chosen from the range 30-150 days. For each observatory, 
we construct a realistic measurement-scatter distribution as 
a function of apparent magnitude (including the blended 
light contribution). For OGLE and MOA, we used represen- 
tative events from the OGLE-III 2009 season (the last season 
of OGLE III), and from the MOA-II 2010 season (obtained 



from the MOA databasqj) , respectively. The scatter was ob- 
tained from the differences between the observed photome- 
try and the best point-source-point-lens fit to each event. For 
Wise and PTF, we used the MOA-II scatter which, based on 
the pilot results, is similar for all three observatories. This 
procedure naturally accounts for all the uncertainties that 
arise from the observations and the data reductions, includ- 
ing outlier photometric measurements (usually due to DIA 
residuals) . Figure [3] (right panels) shows four examples of 
final, sampled and noised, light curves. 



3.5 Anomaly detection 

The most general solution of a microlensing event that has 
an anomaly caused by just one planet requires fitting for 
over a dozen parameters, or even more in the case of multiple 
planets. These parameters represent the source-lens config- 
uration, the planetary system, the dynamics of the system 
(i.e, parallax, orbital motion, xallarap), the source size, and 
the limb-darkening profile. There are degeneracies between 
some of these parameters that can sometimes cause poor 
characterization of the system or even the false detection of 
a planet. 

However, in order to estimate the abundance of plan- 
etary systems, we do not need to fully characterise every 
event, but only to find the signature of a planetary pertur- 
bation with respect to a point-lens light curve. Since the 
mass ratio between a planet and its host star, qt, is rela- 
tively small, the deviation will generally be localized in time, 
when the source passes close to, or crosses, caustics (cen- 
tral or planetary), while the rest of the light curve will stay 
mostly unperturbed (more precisely, the time scale of the 
full light curve will change with respect to a planet-less sys- 
tem, especially for higher-mass ratios, but the form of most 
of the light curve remains that of a point lens) . After fitting 
a global point-lens model, we therefore search each simu- 
lated light curve for a localized deviation from that model, 
using a local x^ test. 

Since a x^ test requires reliable error estimates, we first 
obtain an empirical estimate of the photometric errors, from 
the scatter of the measurements around the point-lens model 
for each light curve of every observatory's dataset, as a 
function of event magnitude. However, when we previously 
noised our simulated light curves (Section 13.41 above) , we 
included also the outlier population of measurement errors 
that exists in the real data. Using an iterative o"-clipping 
procedure (i.e., we repeatedly fit a Gaussian to the error 
distribution, and reject 2a outliers) , we now eliminate these 
outliers and find the best-fitting Gaussian error distribution 
(irrespective of the reported errors, which are often under- 
estimated). With our derived errors, the x^ of the entire 
light curve will naturally be close to the number of degrees 
of freedom. However, when applied piecewise to sections of 
the light curve, x^ can reveal the presence of a planet via 
the temporal condensations of 'outliers' that the planetary 
deviation causes. We define a "running" local x^ estimator: 



Alocal 



= E 



{.h - Uif 



(4) 
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Figure 3. Examples of simulations with our light-curve generator. Left column: Several configurations of the lens plane, in units of 
the lens star's Einstein radius, of projected Solar-like planetary systems around a lens star with a snowline power-law index of s = 1. 
The projected circular orbits of some of the planets are shown (black ellipses). The red circle is the host-star Einstein radius and the 
symbols mark the projected positions of the 8 planets of the system, some of them outside the area shown (Mercury - blue point, Venus 
- blue circle, Earth - brown triangle. Mars - pink diamond, Jupiter - gray square, Uranus - light-gray triangle). The green line is the 
trajectory of the source star and the black cross is the host star. Middle column: Light curves (blue) obtained by ray tracing from the 
observer, through the planetary system, to the source. Black curves are for the lens star without planets. Right column: Light curve for 
each system, as 'observed' by the four observatories of the generation-II microlensing network (green - MOA, blue - OGLE, red - Wise, 
magenta - PTF). The black curves are the light curves for the lens star without any planets. 



where fi and Ui are the "observed" flux and the derived 
error of each epoch, respectively, and /pi is the point-lens 
model flux at that time. The local x^ is repeatedly calcu- 
lated by advancing one observed epoch at a time. We em- 
pirically found the 31-consecutive point length to be opti- 
mal for detecting short-period deviations without increasing 
false-positive detection due to random error correlations. 
An anomaly is considered detected if the y^ probability is 
smaller than 0.3% {x^ > 53.3 for 28 degrees of freedom - 
31 points minus three fitted parameters: ts, uo and to). In 
a real experiment, the method provides also a second-order 
calibration between observatories, if the global photometric 
scatter (for the entire light curve) of a specific observatory 
is not symmetric around zero. 

Photometric outlier points that survived the previous 
filtering will affect the local x^ test and lead to false de- 
tections. In order to avoid this, our anomaly detection cri- 
terion also requires at least three consecutive data points, 
each with a 3a deviation from the point-lens model (as in 



the OGLE E arly-Early Warning System anomaly detector, 
IUdalski|[2003l ). 

Naturally, a more sophisticated or flexible detection cri- 
terion could be devised, one that adapts to different ca- 
dences, weather gaps, and deviation timescales. However, 
for our present purposes, the main point is that, as long as 
the same detection criteria are used on the real and the sim- 
ulated data, the statistical properties of planetary systems 
can still be correctly estimated. 



3.6 Monte Carlo simulation 

For every combination of physical and observational param- 
eters that we study, we generate 3333 simulated light curves, 
as described above, pass them through our detection filter, 
and note the number of planet detections. Naturally, if only 
a fraction / of lens stars host planetary systems, the num- 
ber of detections will be reduced correspondingly, by a factor 
/. We examine combinations of three input parameters in 
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each set of simulation: s - the snowhne scahng power-law 
index; robs - the total observing period, from 1 month (sim- 
ilar to the generation-II pilot campaign of 2010) and up to 
150 days, which is approximately the total visibility period 
of the Galactic bulge every year; and the minimal cadence, 
from 15 minutes, as in the high-cadence fields of OGLE and 
MO A, to 3 hours. 

We recall that the minimal cadence is just the base- 
line for the real observational sequence of the simulation, 
which uses the actual generation-II pilot campaign observ- 
ing times. Similarly, the total observing period includes the 
real gaps (due to, e.g. weather, full moon, technical down- 
time) of the observations. 

The number of expected detections in a generation-II 
seasonal campaign is simply the number of events times the 
detection fraction, and times the true planet abundance, /. 
The number of microlensing events per season in the Galac- 
tic bulge detected by OGLE-III and MOA-II has been about 
600-650. With OGLE-IV, the event rate wiU exceed 1000 
events per season. The fraction of ongoing events inside the 
8 deg^ footprint that was continuously monitored by the 
generation-II pilot has been 34%, so, in what follows, the 
number of expected planetary detections assumes 340 events 
per season. The actual number of detections in a real ex- 
periment will be distributed according to Poisson statistics 
around the expectation values that we predict. 



3.7 Binary contamination 

A potentially important contaminant of a microlensing 
planet survey is "false positives" due to binary, as opposed 
to planetary, systems. With well-sampled, high S/N data, 
binary and planetary deviations can generally be distin- 
guished because of the different secondary/primary mass ra- 
tios. As already noted, the ratio between the timescales of 
the anomaly and the full event is of order the square root of 
the mass ratio. However, in a real generation-II experiment, 
the binary or planetary nature of some events may remain 
ambiguous. To quantify this source of systematic error, we 
simulate a population of binary lenses. We note that, since 
the mass ratio is no longer very small in the case of binaries, 
the perturbation due to the companion mass is no longer lo- 
cal and the single-lens model that best fits the data, before 
applying the detection criteria, will generally be different 
from that of a single lens with the mass of the primary. 

For the primary-mass stars in the simulated binaries, 
we use OGLE-III events as before, and the same mass, 
velocity and distance distributions. We assume that the 

J ihysical binary sep aration, a, is distributed as /(a) ~ a~^ 
Opik fc Lukkl 19241) although o ther forms can be considered 



fe.g. iDuguennov fc MavoiJll99ll ). The mass ratio q of the bi- 



nary is drawn from a flat distribution, f{q) — const, which 
appears to be ge neric to a wide range of binary masse s and 
separations fe.g. ISoderhjeimll2000l . iMazeh et aLlbOOSD . The 
binary contamination in a microlensing survey will naturally 
depend both on the fraction of stars that are in binari es and 
on th e fraction of stars that host planetary systems. iLadal 
l|2006r ) has summarized several studies on binarity fraction. 
He emphasizes that stellar multiplicity is a rising function 
of the stellar mass, and therefore most stars are single (since 
low- mass stars are the most common ones). For example, the 
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Figure 4. Event distribution as a function of total duration of 
a 3(T deviation from a point-lens model, in units of the Einstein- 
radius crossing time, for binaries (black) and for planets (red). 
Planetary deviations are typically much shorter than those due 
to binaries. 



single-star fraction is ~70% among M-type stars and ~45% 
among G-type stars. 

Figure |4] shows, for the simulated populations of plan- 
ets and binaries, the distributions of total detected devia- 
tion time (i.e. the total time there was a 3a deviation from 
the best-fit point-lens light curve). The deviation time is in 
units of ts for each event. As expected, the typical binary 
deviations are of order ts while the most common planetary 
deviations are two orders of magnitude more brief. However, 
there is an overlap in the distributions, which can be a source 
of false positive planet detections. For example , assuming a 
planet ary frequ ency o f / = 1 /6, as estimated bv lGould et al.l 
(|2010f ). and the iLadal 1 20061 ) binarity fraction as a function 
of primary mass, the contamination of binaries to the short 
deviation events (total deviation less than O.Ue) would be 
33%. 

The binary population is currently much better char- 
acterised than the planetary one. Since binaries completely 
dominate the events with perturbations having timescales 
> O.Ue, that part of the distribution can be used to anchor 
a model for the binary population, and thus to estimate and 
statistically subtract off the binary contamination on short 
timescales. Furthermore, although our emphasis has been on 
the detection of deviations from single-lens light curves, de- 
tailed modeling that discriminates between planets and bi- 
naries will certainly be possible for a fraction of the events. 
Such modeling will be possible either with the Generation-II 
data themselves, or with the addition of foUowup data from 
networks that focu s on anomalous e vents, such as Micro FUN 
(|Gould et al.ll201Gl ) and PLANET l|Albrow et al.lll998l ). 



4 RESULTS 

In this section, we present the results of our simulations for 
various observational and physical input parameters. 
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Figure 5. Le/t: Number of planetary detections as a function of the minimal cadence for different snowline indices s, for a total observing 
period of 150 days, with 340 events per season in the generation-II survey footprint, and assuming each lens star hosts a solar-analog 
planetary system. Multiple-planet detections are counted as one event. The right-hand vertical axis is the detection fraction among all 
of the events monitored. For a fraction / of lens stars hosting planets, the curves will scale down correspondingly. Right: Planetary 
detections as a function of the total observing time per season, T^hsy with a minimal cadence of 15 minutes for different snowline indices 
s, as before, assuming each lens star hosts a planetary system. 
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Figure 6. Left: Number of detections as a function of peak /-band magnitude of an event, for a 150-day season and various cadences, 
with results averaged over all snowline indices. The maximum is at ^17 mag. Right: Detected fraction of events versus peak /-band 
magnitude. The increase towards brighter events results both from the more accurate photometry and from the higher sensitivity to 
planets of high-magnification events. 



4.1 Observational parameters 

Figure [5] shows the total number of planetary detections 
per season as a function of the minimal observing cadence, 
assuming every lens star hosts a planetary system, a to- 
tal observing period of 150 days and for a range of snowline 
power-law indices, s. We find that over half of the detections 
are missed by going from a minimal cadence of 15 minutes to 
3 hours. This is simply a consequence of the loss of low-mass 
planet detections, with their short timescales. Figure [5] also 
shows the total number of planetary detections as a func- 
tion of the total observing period, assuming each lens star 
hosts a planetary system, a minimal cadence of 15 minutes 



and the various snowline indices. As expected, the number 
of detections rises, more or less linearly with the total ob- 
serving time, because the total number of events is propor- 
tional to the total observing period. In addition, there is a 
second-order contribution to the detected number for long 
experiment durations because a detectable deviation in long 
events (long ts) can be missed by a brief campaign. 

The detection fraction is naturally sensitive to the pre- 
cision of the photometry, which depends on the event magni- 
tude. As mentioned in Section [ij this is already evidenced in 
the planets that have been discovered to date via microlens- 
ing, which have typical peak event magnitudes of ~15. These 
events were either of high-magnification, or with sources that 
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Figure 7. Planet detection efficiency as a function of impact 
parameter, for a 150-day season and various cadences, with results 
averaged over all snowline indices. 



were bright in the first place, and with long event time scales, 
tE- Generation-II experiments improve upon this limitation 
with high cadences and accurate photometry. Figure|6]shows 
that the number of planetary detections as a function of I- 
band magnitude at peak, for the type of events that are 
detected by the current surveys (OGLE and MO A), has 
a maximum at magnitude ~17 for a 150-day season and 
various cadences, averaged over all snowline indices. The 
peak results from, on the one hand, the increased number of 
events with low magnifications and faint baselines, and on 
the other hand, the reduced sensitivity of low-magnification 
and faint events. The detected fraction of events per peak 
7-band magnitude increases about linearly with magnitude 
for brighter events (Figure[6l right panel). This folds in both 
the high photometric precision of the bright events and the 
higher sensitivity to planets of high-magnification events, 
which we further explore in the next section. 



4.2 Physical parameters 

The sensitivity to planets in high-magnification events is 
close to 100% for planets projected near the Einstein ra- 
dius, because of their s ignatures on the central caustic 
l|Griest fc Safizadehll99a ). The sensitivity naturally declines 
for lower impact parameters. Our simulations permit us 
to study this question for the specifics of the generation- 
II experiment. Figure [7| shows the detection fraction of 
events as a function of the impact parameter uo (with re- 
sults averaged over all snowline indices). The detected frac- 
tion depends on impact parameter approximately as Uq^ . 
Among intermediate magnification events, with impact pa- 
rameters 0.1 < Uo < 0.3 (corresponding to amplifications 
10 > Aniax > 3.4) planetary systems are detected with 
~17% efficiency for the optimal observing setup (150-day 
season and minimal cadence of 15 minutes). 

Figure |8] shows the number of planetary detections as 
a function of s, for a total observing season of 150 days, 
and for various cadences. The detection is insensitive to the 
snowline index in the range 0.5 ^ s ^ 1, but there is a drop 
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Figure 8. Planetary detections as a function of the snowline 
scaling index s, -Rsnow oc M", for various cadences. As in Figs.[5]- 
\E\ all stars host solar-analog systems, and a 150-day bulge season 
is assumed. 



at s > 1 such that, by s = 2, the number of detections is 
lower by about 40% compared to s = 0.5. Low values of s 
keep the 'Jupiters' in our lens systems close to the Einstein 
radius, which scales with mass to the 1/2 power. Conversely, 
larger values of s shift the giant planets away from the Ein- 
stein radius, where they are less likely to produce detectable 
anomalies (see Figure |9]l . The rocky planets that are then 
sometimes near the Einstein radius have masses that are too 
low for a probable detection. 

As a concrete ex ample, for th e frequ ency of Solar-like 
systems estimated bv lGould et al.l (|20ld ). of / = 1/6, and 
assuming a snowline index of 1, the number of expected 
planet detections per season is 11.3 ± 3.4 in an optimal 
generation-II experiment. A four-year generation-II experi- 
ment can thus potentially produce a sizeable sample of order 
50 microlensing detected planets. Conversely, for an assumed 
value of s, such an experiment will determine the abundance 
of snowline planets to a precision of ~ 15%, or even better if 
the actual planet frequency is higher. Naturally, this result 
folds in our assumed scaled solar-analog model, with its par- 
ticular range of planetary locations. This assumption leads 
to a predominance of gas and ice giants near the Einstein 
radius, though all types of planets can be near the Einstein 
radius due to projection elfects (see Figure [9}. 



4.3 Multiple planet detection and lensing planet 
identity 

For the optimal observing plan of 150 days with 15-minute 
minimum cadence, we find that 83% of the detections are 
due to deviations produced by a single planet, 16% are due 
to two planets (i.e., two planets contribute to the detected 
anomalies in an event) and fewer than 1% of the detections 
are caused by three of the planets in the system. Thus, if 
solar analogs are common, one out of every 4 or 5 planetary 
detections could tu rn out to be a double-planet detection, a s 
in OGLE-06-109L (|Gaudi et al.ll2008l : iBennett et al.ll2010l ). 
In terms of the identity of the planets causing the anoma- 
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Figure 9. Left: Mean projected planetary separations of the simulated Jovian planets in units of the Einstein radius, r^ = 6e ■ Dl, 
as a function of the snowline scaling index. Error bars indicates the RMS of each distribution. The gray area is the "lensing zone" , 
0.6 ^ t/te ^ 1.6, where microlensing is most sensitive to planets. Right: Number of detections of 'Jupiters' , 'Saturns' and 'Neptunes' 
as a function of the snowline scaling index, for the optimal observing plan of 150-day season and minimal cadence of 15 minutes. 



lies, for a scaling index s — 2, over 65% of the planets are 
'Jupiters', 24% are 'Saturns', 4% 'Uranus', 4% 'Neptunes', 
and 1% 'Earths'. The fraction of detected 'Jupiters' rises 
for lower snowline scaling indices, reaching 84% for s — 0.5. 
Figure [9] (left panel) shows the projected locations of the 
Jovian planets in units of the Einstein radius, ve = Oe ■ Dl, 
as a function of the snowline scaling index. Microlensing 
is most sensitive to planets in the so-calle d "lensing zone", 
0.6 < r/rE ^ 1.6 (JGriest fc Safizadch"l998''). The right panel 
of Figure [9] shows the number of detections of Jupiters, Sat- 
urns and Neptunes as a function of the snowline scaling in- 
dex. As seen in this figure, the lensing dominance of Jupiters 
among the detections is a consequence of the fact that, un- 
der our scaled solar-analog assumption, this planet (the most 
massive one) is never far from the Einstein radius. Never- 
theless, medium-mass planets, such as Neptune, can also be 
detected with high cadence surveys, whenever they are pro- 
jected near the lensing zone. 



5 SUMMARY 

In this work, we have studied various aspects of generation-II 
microlensing experiments - the optimization of the observa- 
tional parameters of such experiments, and their sensitivity 
to planets for different physical assumptions about the plan- 
etary systems. Our main results are the following: 

(i) Long duration experiments with continuous high ca- 
dence are vital. Experiments with low cadences of 3 hours 
will detect only about half as many planets as experiments 
with cadences of 15 minutes, and likewise experiments last- 
ing one-half, versus a full, bulge season. 

(ii) The current approach of using a network of 1-m-class 
telescopes for a generation-II experiment is adequate for 
tracking the low-magnification events discovered by these 
same telescopes. We have quantified, for our particular net- 
work and its parameters, the decrease of the planet detec- 



tion efficiency with the decrease of the peak brightness of an 
event. 

(iii) The abundance of planetary systems is degenerate 
with the snowline scaling index. Conversely, with constraints 
on planetary fraction from other planet-hunting techniques, 
microlensing surveys could discriminate among different 
planet formation scenarios that predict different snowline 
scalings. 

(iv) The majority of events with anomalies will be due 
to single planet deviations (83%), assuming Solar system 
distance and mass ratios, but the non-negligible remainder 
can reveal two planets. 

(v) For Solar-like systems, at least 65% of the detections 
will be of 'Jupiters', due to both the mass and the location 
near the Einstein radius of this planet. 

(vi) The current generation-II experiment which continu- 
ously monitors 8 deg^ of the Galactic bulge with the highest 
microlensing event rate, and with a cadence of 15-30 min- 
utes, will be able to determine, after four full seasons, the 
abundance of planets near the snowline to a precision of 10- 
30%. This, assuming the fraction of Galactic stars hosting 
such planets is between 1/3 to 1/10, and a snowline scaling 
power-law index in the range 0.5 to 2. 

(vii) The main contaminant to planet detection, bina- 
ries, can largely be filtered out already at the detection 
stage, based on the longer deviation timescales that they 
produce. Additional characterization and modeling of indi- 
vidual events in the timescale overlap region between planets 
and binaries will further lower the contamination fraction. 

A full generation-II experiment, of the type we have 
simulated here, is already in progress, our results show that 
first estimates of the frequency and nature of snowline- region 
planets will be possible soon. To improve the analysis, our 
future simulations will include also higher-order dynami- 
cal effects (parallax, orbital motion, xallarap), a distribu- 
tions of source sizes, limb darkening, and a wider range of 
physical planetary system possibilities. Together with re- 
sults from other planet-search techniques, a clearer picture 
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should emerge of the properties of planetary systems and 
the physics of their formation. 
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